arXiv: 1503.08243vl [cond-mat.quant-gas] 27 Mar 2015 


Observation of Bose-Einstein Condensation in a Strong Synthetic Magnetic Field 


Colin J. Kennedy*, William Cody Burton, Woo Chang Chung, and Wolfgang Ketterle 

MIT-Harvard Center for Ultraeold Atoms, Researeh Laboratory of Eleetronies, 
Department of Physies, Massaehusetts Institute of Teehnology 
(Dated: March 31, 2015) 

Extensions of Berry’s phase and the quantum Hall effect have led to the discovery of new states of 
matter with topological properties. Traditionally, this has been achieved using gauge fields created 
by magnetic fields or spin orbit interactions which couple only to charged particles. For neutral 
ultracold atoms, synthetic magnetic fields have been created which are strong enough to realize the 
Harper-Hofstadter model. Despite many proposals and major experimental efforts, so far it has not 
been possible to prepare the ground state of this system. Here we report the observation of Bose- 
Einstein condensation for the Harper-Hofstadter Hamiltonian with one-half flux quantum per lattice 
unit cell. The diffraction pattern of the superfluid state directly shows the momentum distribution 
on the wavefuction, which is gauge-dependent. It reveals both the reduced symmetry of the vector 
potential and the twofold degeneracy of the ground state. We explore an adiabatic many-body state 
preparation protocol via the Mott insulating phase and observe the superfluid ground state in a 
three-dimensional lattice with strong interactions. 


Topological states of matter are an active new frontier 
in physics. Topological properties at the single particle 
level are well understood; however, there are many open 
questions when strong interactions and correlations are 
introduced mm as in the = 5/2 state of the frac¬ 
tional quantum Hall effect [3] and in Major ana fermions 
[M]. For neutral ultracold atoms, new methods have 
been developed to create synthetic gauge fields. Forces 
analogous to the Lorentz force on electons are engineered 
either through the Coriolis force in rotating systems 13 - 
[5], by phase imprinting via photon recoil [TUHIl], or 
lattice shaking HU [IS]. Much of the research with ul¬ 
tracold atoms has focused on the paradigmatic Harper- 
Hofstadter (HH) Hamiltonian, which describes particles 
in a crystal lattice subject to a strong homogeneous mag¬ 
netic field [HHini- For magnetic fluxes on the order of one 
flux quantum per lattice unit cell, the radius of the small¬ 
est possible cyclotron orbit and the lattice constant are 
comparable, and their competition gives rise to the cel¬ 
ebrated fractal spectrum of Hofstadfer’s butterfly whose 
sub-bands have non-zero Chern numbers m and Dirac 
points. 

So far, it has not been possible to observe the ground 
state of the HH Hamiltonian, which for bosonic atoms 
is a superfluid Bose-Einstein condensate. It is unknown 
whether this is due to heating associated with technical 
noise, non-adiabatic state preparation, or inelastic colli¬ 
sions. These issues are complicated, since all schemes for 
realizing the HH Hamiltonian use some form of tempo¬ 
ral lattice modulation and therefore are described by a 
time-dependent Floquet formalism. The HH model arises 
after time averaging the Floquet Hamiltonian, but it is 
an open question to what extent finite interactions and 
micromotion lead to transitions between Floquet modes 
and therefore heating [2T]-[24]. Bose-Einstein condensa¬ 
tion has been achieved in staggered flux configurations 
[I5l|25] and in small ladder systems |27] further high¬ 


lighting its noted absence in the uniform field configura¬ 
tion. 

In this article, we report the first observation of Bose- 
Einstein condensation in Hofstadter’s butterfly. The 
presence of a superfiuid state in the HH lattice allows 
us to analyze the symmetry of the periodic wavefunction 
by self-diffraction of coherent matter waves during ballis¬ 
tic expansion - analogous to a von Lane x-ray diffraction 
pattern, which reveals the symmetry of a lattice. Using 
this method, we show that the unit cell of the magnetic 
lattice is larger than that of the underlying cubic lattice. 
This reflects that the vector potential necessarily has a 
lower symmetry than a uniform magnetic field, since a 
transitionally invariant field can only be realized by a vec¬ 
tor potential that breaks this translational invariance. In 
our experiment, we analyze the wavefunction directly and 
observe both its non-gauge invariance [28] and ground 
state degeneracy. Finally, we explore a many-body adi¬ 
abatic state preparation protocol where the many-body 
gap of the Mott insulator is used to switch superfiuid 
order parameters without creating excess entropy. We 
realize the HH Hamiltonian superfiuid ground state in a 
three-dimensional lattice with strong interactions. 

The HH model with large synthetic magnetic fields is 
realized as an effective Hamiltonian engineered by laser- 
assisted tunneling processes in a tilted lattice potential. 
As in Refs. dSKH], tunneling in the x-direction of a two- 
dimensional optical lattice is suppressed by an energy 
offset and then subsequently restored with a resonant 
Raman process with a momentum transfer, 6k = kxX -\- 
kyp. The z-direction consists of loosely confined tubes 
of condensate, so interactions are weak. In the time- 
averaged picture, this gives rise to the Hamiltonian: 

H = (^ — Ke ~ J 

( 1 ) 

where (frn^n = Tnkxa + nkya = 7r{m + n) reflects the spe- 
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FIG. 1. Observation of Bose-Einstein condensation in the Harper-Hofstadter model with one half flux quantum per unit cell, 
(a) Spatial structure of the cubic lattice with the synthetic vector potential - (dashed) x-bonds feature a spatially dependent 
tunneling phase 0mn = 7r(m + n), whereas tunneling along (solid) y-links is the normal tunneling. The synthetic magnetic 
field generates a lattice unit cell that is twice as large as the bare cubic lattice (green diamond), (b) The band structure 
of the lowest band of the HH Hamiltonian shows a twofold degeneracy of the ground state and two Dirac points located at 
ky — ±7r/2a. The magnetic Brillouin zone (green diamond) has half the area of the original Brillouin zone. Due to degeneracies 
related to the indistinguishability of the sites in the magnetic unit cell, the primitive cell of the band structure is even smaller 
(doubly reduced Brillouin zone, brown square). These lattice symmetries are both revealed in time-of-flight pictures showing 
the momentum distribution of the wavefunction. (c-f) Schematics of the momentum peaks of a superfluid. The dominant 
momentum peak (filled circle) is equal to the quasimomentum of the ground state. Due to the spatial periodicity of the 
wavefunction, additional momentum peaks (empty circles) appear, separated by reciprocal lattice vectors (green arrows) or 
vectors connecting degenerate states in the band structure (brown arrows), (g-j) Time-of-flight images. The superfluid ground 
state of the normal cubic lattice is shown in (g) in comparison to the superfluid ground state of the HH lattice (h-j). In (h), 
only one minimum of the band structure is filled, directly demonstrating the fundamental symmetry of the HH Hamiltonian 
in our chosen gauge. The number of momentum components in (i-j) is doubled again due to population of both degenerate 
ground states of the HH Hamiltonian. The micromotion of the Floquet Hamiltonian is illustrated in (e-f, i-j) as a periodically 
shifted pattern in the x-direction, analogous to a Bloch oscillation. All diffraction images were taken after at least 30 ms hold 
in the HH lattice. The field of view of all experimental pictures is 631 /im x 631 /xm. 


cific gauge implemented (Fig. la). An atom that trav¬ 
els around a single lattice unit cell picks up a phase of 
TT = kya^ with a being the lattice spacing, so the Hamil¬ 
tonian is equivalent to that of a charged particle in a 
magnetic field with one half of a magnetic flux quantum 
per unit cell. 

After preparing the atoms in the HH lattice (see Meth¬ 
ods), we observe the momentum distribution of the wave- 
function by suddenly turning off all laser beams to allow 
20 ms time-of-flight, and measure the resulting density 
distribution with absorption imaging (Fig. Ih-j). We 
first note that the images show sharp peaks, which are 
the hallmark of a superfluid Bose-Einstein condensate in 
a periodic potential. This demonstrates for the first time 
that we have successfully prepared a low entropy state 
in the bulk HH Hamiltonian. In addition, the diffrac¬ 
tion images directly show a reduced symmetry of the su¬ 
perfluid, despite the translational symmetry of both the 
lattice and the homogeneous synthetic magnetic field. 


Fundamentally, the vector potential is responsible for 
the broken translation symmetry of the lattice, and the 
time-of-flight patterns can be understood by examining 
these symmetries. In real space, there are two relevant 
unit cells: one is the unit cell of the cubic lattice, and 
the other is the unit cell of the Hamiltonian determined 
by the experimental gauge, which we call the magnetic 
unit cell. In highly symmetric gauges - such as our ex¬ 
perimental gauge and the Landau gauge - a magnetic 
flux oi a = pjq has a magnetic unit cell that is q times 
larger than the original unit cell and contains q indistin¬ 
guishable sites [28]. In momentum space, these unit cells 
correspond to the Brillouin zone of the underlying lattice 
and to the magnetic Brillouin zone, respectively. 

Due to the indistinguishability of the sites in the 
magnetic unit cell, every quasimomentum state within 
the magnetic Brillouin zone is g^-fold degenerate with 
other quasimomentum states connected by the transla¬ 
tion symmetry operators of the original lattice modified 
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by a phase 27r/g. This modified translation symmetry 
connects all degenerate ground states in the magnetic 
Brillouin zone, and corresponds to a third relevant scale 
2711 (qa) useful for understanding the diffraction pattern 
of a superfluid in the HH Hamiltonian. This length scale 
defines the doubly reduced Brillouin zone, which is the 
primitive cell of the band structure [28]. 

Experimentally, we observe the gauge-dependent 
diffraction pattern directly when the condensate occu¬ 
pies only one minimum of the band structure (Fig. Ih). 
The pattern’s symmetry reflects that of the gauge po¬ 
tential (Fig. Id), directly demonstrating the symmetry 
of the reciprocal lattice vectors a{x±y) of the magnetic 
Brillouin zone. A common belief is that all observables 
are gauge-independent. However, gauge-dependent ob¬ 
servations can be made in time-of-flight images of ul¬ 
tracold atoms when the momentum distribution of the 
wavefunction is observed. The sudden switch-off of all 
laser beams transforms canonical momentum, which is 
gauge dependent, into mechanical momentum, which is 
readily observed |29| . 

The diffraction pattern shown in Fig. li reveals the 
symmetry of the band structure (Fig. lb). In the exper¬ 
imental image, both degenerate states are equally popu¬ 
lated and the peaks can be regarded as the sum of the 
diffraction pattern in Fig. Ih with an additional copy 
separated by the gauge-invariant scale 27r/(ga) = 7r/a in 
either x or y directions. The diffraction pattern reflects 
this two-fold multiplicity of momentum peaks in both the 
X and y directions and features a four-fold increase of vis¬ 
ible peaks m- To the best of our knowledge, this work 
is the first direct observation of the different symmetries 
characterizing a simple lattice with a homogeneous mag¬ 
netic field: Figs. Ig-i show the symmetry of the lattice, 
the gauge potential, and the band structure. Realization 
of a Bose condensed ground state allows direct obser¬ 
vation of both gauge-dependent and gauge-independent 
symmetry breaking. 

A final feature of the observed diffraction pattern is the 
periodic oscillation of the peaks along the x-axis. This 
is the result of micromotion coming from acceleration 
due to the tilt and Raman drive. In the limit of a low 
drive amplitude, this micromotion is a Bloch oscillation 
in the tilted lattice. More generally, it is the effect of 
a time-dependent transformation between the rotating 
frame in which the HH Hamiltonian is realized and the 
lab frame in which it is observed mi. Figs. Ie,i and 
f,j show diffraction patterns for the micromotion at zero 
phase and at 7r/2 phase - equivalent to zero and tt phases 
in a Bloch oscillation due to the doubling of the unit cell. 
In the experiment, the absolute phase of the micromotion 
is a random variable due to wavelength scale drifts of the 
relative positions of the stationary and Raman lattices. 

Qualitatively, almost all images show roughly equal 
population of both degenerate ground states. Occupa¬ 
tion of only one state (as in Fig. Ih) is observed in less 




FIG. 2. Population imbalance of the two ground states of the 
HH Hamiltonian with 1/2 flux, (a) Band mapping sequence 
adiabatically connecting quasimomentum to free space mo¬ 
mentum. The Raman beams were ramped down in 0.88 ms, 
followed by a linear ramp down of the lattice beams in 0.43 ms. 
(b) The histogram shows the relative population imbalance 
of the two degenerate minima, consistent with a spontaneous 
symmetry breaking mechanism during state preparation. The 
data consists of 30 shots taken after a hold time of 29.4 ms 
in the HH lattice. The inset displays a raw image for the 
band mapping of the 1/2 flux superfluid with two degenerate 
ground states compared to a topologically trivial superfluid 
with one ground state. 

than 1% of the shots. For a quantitative determination 
of the populations, we developed a band mapping tech¬ 
nique that adiabatically connects the HH ground states 
to momentum states (see Supplement). Band mapping 
reduces the complicated diffraction patterns of Fig. I to 
two peaks, one for each degenerate ground state. Fig. 
2b shows that the minima are most frequently populated 
with equal proportion, demonstrating the degeneracy of 
the minima and the robustness of the loading procedure 
to technical fluctuations. However, mean-field calcula¬ 
tions predict that a superposition state, which would 
create spatial density modulation, is energetically unfa¬ 
vorable due to repulsive interactions [28]. Therefore, the 
simultaneous occupation of both minima most likely in¬ 
dicates the presence of domains, composed of atoms in 
only one of the dispersion minima. 

To rule out non-adiabatic transitions to higher bands 
as an alternative explanation of the two peaks in band 
mapping, we tested the technique on a topologically triv¬ 
ial Floquet superfluid (see Methods). This system con¬ 
sists of a cubic lattice with a tilted potential where tun¬ 
neling in the tilted axis is restored by amplitude modula¬ 
tion (AM) [311 [32]. Its micromotion, which is similar to 
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FIG. 3. Decay of Bose-Einstein condensates in modulated 
lattices. The figure compares the decay of the 1/2 flux HH 
superfluid (red circles) to the decay of the AM superfluid (blue 
squares). Note that the lower visibility of the HH superfluid is 
due to the peak doubling, which at the same condensate frac¬ 
tion, leads to lower visibility. Exponential fits to the decay of 
the visibility of the diffraction patterns give lifetimes of 142 ± 
15 ms and 71 ± 8 ms, respectively. The data were taken after 
a 10 ms hold time after switching on the final Hamiltonian us¬ 
ing the non-adiabatic procedure (see Methods). Uncertainty 
is given by the statistical error of the mean of five repetitions 
of the experiment, added in quadrature to uncertainty in the 
peak visibility fitting. 


that of the HH system, and its trivial band structure al¬ 
lows identification of non-adiabatic excitations attributed 
to the micromotion. The resulting momentum distribu¬ 
tions for both the AM and the HH systems are shown 
in Figure 2b. Non-adiabatic band excitations would ap¬ 
pear as additional peaks in the tilted direction, which are 
not observed. Therefore, we interpret the band mapping 
procedure in the HH system as an adiabatic process. 

The above results confirm many aspects of the weakly 
interacting superfluid ground state of the HH Hamilto¬ 
nian. However, there are many open questions about 
topological physics in the presence of strong interactions, 
often involving small energy differences. To study this 
regime, long coherence lifetimes are necessary. 

In order to evaluate our coherence lifetime and disen¬ 
tangle different sources of decoherence and heating, we 
extensively used the AM superfluid as a benchmark. Ini¬ 
tially, the HH superfluid had a vastly shorter lifetime 
than the AM superfluid. Subsequent improvements in 
the stability of magnetic fields, beam pointing, Raman 
phase, and gradient alignment to the lattice direction 
improved the HH lifetime until both systems had a com¬ 
parable lifetime in the ground state, given in Fig. 3. The 
lifetimes for the number of trapped atoms is much longer 
than the coherence lifetime in either case, with no dis¬ 
cernible loss of atom number up to 500 ms. This leads 
to the conclusion that two-body dipolar interactions and 
three-body recombination, which would lead to particle 


loss from the trap, are not limiting the lifetime of the 
ground state. The same applies to excitations to higher 
lying bands which, due to the tilt, are strongly coupled 
to the continuum via Landau-Zener tunneling. Therefore 
the decay in Fig. 3 is dominated by transitions within the 
lowest band, which can be caused by technical noise or 
by elastic collisions that transfer energy from the micro¬ 
motion into heat. While our heating seems to be largely 
technical, it is useful to briefly discuss these collisions as 
an ultimate limit on the lifetime. 

For the tilted lattice, one decay path is via overlap 
of neighboring Wannier-Stark states with offset energy 
A, which can be transferred to excitations of the low¬ 
est bands, or to the free particle motion along the tubes 
orthogonal to the 2D lattice. Such processes are de¬ 
scribed by Fermi’s golden rule for transitions between 
Floquet states Emu- Using this framework, we derive 
the scattering rate for the AM superfluid to be between 
F ^ 0.30—0.68 s“^, where the uncertainty comes from the 
uncertainty in our density measurement due to redistri¬ 
bution during lattice ramp-up. This estimate is smaller 
than the observed decay rate by a factor of ^ 10 — 20. 
See Supplement for a derivation and experimental pa¬ 
rameters. 

In principle, an adiabatic procedure to prepare the 
ground state of the HH Hamiltonian should enable bet¬ 
ter control and higher fidelity in the final state. How¬ 
ever, we empirically found that a sudden turn-on of the 
tilt and the Raman beams gave the most consistent high 
contrast images, and was therefore useful in evaluating 
technical improvements. A slower turn-on could intro¬ 
duce more technical heating, and the entropy created in 
the sudden turn-on could be absorbed by other degrees 
of freedom. Such an entropy reservoir is provided in the 
2D lattice schemes by the third dimension - in our case 
tubes of condensates with ^500 atoms. For extensions to 
3D lattices, we implement an adiabatic state preparation 
procedure. 

Adiabatic methods usually require a pathway where 
the ground state is protected by a gap. For single par¬ 
ticle states, this may involve matching the size of the 
unit cell between a topologically nontrivial lattice and 
a trivial lattice PUI [53]. Instead, we explore the use of 
many-body gap of the bosonic Mott insulator state. Since 
the Bose-Hubbard model and the HH model with inter¬ 
actions have the same Mott insulating ground state for 
large U (assuming the high frequency limit 7/ ^ A), one 
can connect a trivial superfluid in the cubic lattice to the 
ground state of the HH model via two quantum phase 
transitions. The first transition “freezes” out the origi¬ 
nal superfluid by entering the Mott insulator where the 
phases of the tunneling matrix elements can be changed 
to a non-zero flux configuration without adding entropy 
|33| . The second transition “unfreezes” the Mott state 
into the HH superfluid. Technically, this approach can¬ 
not be fully adiabatic due to the exact degeneracy of the 
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FIG. 4. Experimental sequences for two different state prepa¬ 
ration protocols. The non-adiabatic scheme switches sud¬ 
denly from a standard 2D lattice to the HH lattice. The 
adiabatic protocol uses a quantum phase transition to the 
3D Mott insulator as an intermediate step. For this, lattice 
beams are adiabatically ramped up in all three directions to 
20 Er to enter the Mott insulating phase, after which the hy- 
perfine state is flipped with a 0.29 ms RF sweep as done in the 
non-adiabatic scheme. The lattices are then ramped down to 
their final values while the Raman lattice is ramped up in 35 
ms. Lifetimes of both methods are given in the top corner of 
each figure; the lifetime of the adiabatic approach is measured 
from the end of the 35 ms ramp. Errors are statistical and 
given by the exponential fit of the peak visibility. 


ground states in the HH spectrum; a Kibble-Zurek model 
implies that this will lead to the spontaneous formation 
of domains [34]. However, this will not create entropy 
beyond the randomization of domains. 

We implement this scheme to load the ground state of 
the HH Hamiltonian. So far, we have found the adiabatic 
turn-on to be less robust against technical noise, leading 
to a higher shot-to-shot variability; however, the best 
images have identical visibility and a marginally better 
lifetime than that of the non-adiabatic state preparation 
method (Fig. 4). 

We can now use the lattice in the third dimension, re¬ 
quired for the adiabatic loading protocol, to add stronger 
interactions to the physics of gauge fields. To do this, 
we keep the lattice depths and tunneling strengths along 
the X- and y-directions constant, but vary the depth of 
the z-lattice, thereby adjusting the Hubbard interaction 
parameter U from nearly zero to ^450 Hz. With an es¬ 
timated filling factor of ^5, we can access interaction 
energies in excess of x 2 kHz, which is close to the 
superfluid-to-Mott insulator transition in high magnetic 
fields [35] . 

We realize HH superfluids with lattice depths from 0 
to 20 Er and observe at least weak superfluid diffraction 
peaks (see Supplement). Coherence lifetimes for z-lattice 
depths up to 11 Er are presented in Fig. 5. The aver¬ 
age visibility of the diffraction pattern is reduced at high 
z-lattices, but the lifetimes are reduced only by a factor 


Hubbard Parameter U (Hz) 



FIG. 5. HH model with strong interactions. Shown is the 
lifetime of the visibility of the diffraction pattern versus z- 
lattice depth. The top axis shows the Hubbard interaction 
parameter U. All lifetimes are measured from the end of the 
35 ms ramp exiting the Mott insulator. Uncertainty is given 
by the statistical error of the mean of five repetitions of the 
experiment, added in quadrature to uncertainty in the peak 
visibility fitting. 


of ^ 2 — 4 from the lifetime of the weakly interacting 
superfluid at 0 Er. However, we observe much larger 
shot-to-shot variation: some shots will have clearly vis¬ 
ible diffraction peaks, while other attempts show no co¬ 
herence at all. We attribute the variation to technical 
fluctuations, which seem to have an enhanced effect at 
higher 3D lattices. 


In conclusion, we have observed Bose-Einstein conden¬ 
sation in a strong synthetic magnetic field. We have di¬ 
rectly observed how the gauge potential breaks the sym¬ 
metry of the original lattice. To reach the regime of 
strong interactions, we have demonstrated that the HH 
model can be implemented with reasonably long coher¬ 
ence times in the presence of both strong gauge fields and 
strong interactions. Therefore, this work lays the ground¬ 
work for further investigations of topologically nontrivial 
and/or highly frustrated states. Straightforward modifi¬ 
cations of our setup such as the introduction of complex 
tunneling in the z-direction or pseudospins realize Weyl 
points [36|, the quantum spin Hall Hamiltonian, and ul¬ 
tracold atom topological insulators m- Recent tech¬ 
nological advances in controlling and detecting ultracold 
atoms combined with this work offer a promising path 
toward realization of bosonic Laughlin states in optical 
lattices (35] [38f[40] . The ability to realize strong gauge 
fields with strong interactions highlights the capability of 
ultracold atom systems to realize exotic states of matter 
with no known analogues in nature. 
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METHODS 

The experiment begins with a nearly pure BEG con¬ 
fined in a crossed dipole trap in the |1,—1) hyperfine 
state with ^7 x 10^ ^^Rb atoms. We turn on a magnetic 
field gradient to levitate the atoms against gravity and 
simultaneously weaken the dipole traps by lowering the 
power to their final values. Erom here, the non-adiabatic 
and adiabatic sequences differ. 

In the nonadiabatic sequence, the condensate is adi- 
abatically loaded into a two-dimensional optical lattice 
composed of one vertical beam and one horizontal beam 
with lattice constants 532 nm and depth of 11 recoil ener¬ 
gies (Er) in both the x- and y-directions, in the presence 
of a very weak Raman lattice (< 0.1 Er) with relative 
frequency detunings equal to the Bloch oscillation fre¬ 
quency 3,420 Hz. The two beams are beat together on 
a beam cube located as close to the atoms as our opti¬ 
cal setup allows, and the signal is used to stabilize the 
relative phase of the Raman beams. 

Once a phase lock is achieved, we turn on a large tilt 
in the regime J <C A <C Egap by sweeping the frequency 
of an RE field in 0.29 ms to transfer all the atoms to 
the |2, — 2) hyperfine state. Here J is the bare tunneling 
energy in the lowest band, A = 3,420 Hz is the energy 
offset between adjacent sites, and Egap is the energy gap 
between the lowest band and the first excited band. Upon 
completion of the RE sweep, the initial system represents 
an uncoupled stack of 2D BEG’s. 

Eor topologically trivial superfluids, resonant tunneling 
is re-established with amplitude modulation of the ver¬ 
tical (x-axis) lattice with modulation frequency A. Eor 
the 1/2 flux superfluid resonant tunneling is instead re¬ 
established by linearly ramping up both Raman beams in 
0.58 ms to 2U/A = 1.6, where U is the two-photon Rabi 
frequency. After a variable hold time, all laser beams 
are switched off to allow 20 ms time-of-flight, followed by 
absorption imaging. 

In the adiabatic sequence, the condensate enters the 
Mott insulating phase after lattices in all three direc¬ 
tions are adiabatically ramped up to 20 One of the 


Raman beams is then ramped to its final value in 15 
ms, after which the hyperfine state is flipped with a 0.29 
ms RE sweep as explained previously. The lattices are 
then ramped down to their final values while the second 
Raman beam is ramped up in 35.09 ms, thereby adia¬ 
batically connecting the Mott insulator to the 1/2 flux 
superfluid. Again, time-of-flight absorption imaging is 
performed after a variable hold time. 
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SUPPLEMENT 


1/2 Flux Superfluid Diffraction Pattern Fitting 


We quantify the visibility of the 1/2 flux superfluid diffraction pattern above the broad, thermal background by 
fitting one-dimensional line cuts of the column density in both the x- and y-lattice directions. To design a peak 
visibility fitter for the HH superfluid, one needs to account for multiple degrees of freedom associated with the new 
peaks that emerge, in addition to the trivial superfluid peaks. In addition, the uncontrolled phase of the micromotion 
leads to changes in the diffraction pattern, as seen in the paper. This is because the relative optical path length 
between the cubic lattice and the running wave Raman lattice at the position of the atoms is unstabilized. When 
designing a fitter, we must build in the flexibility in handling shot-to-shot changes in the peak locations on the CCD 
camera. 

The fitter works by finding the center of mass of the image and taking the average of several rows of pixels in each 
direction around this center value. The total visibility is determined by least-squares fitting both the x and y line 
cuts to the sum of six gaussians - one broad thermal background and five sharp diffraction peaks. The visibility of a 
given peak is defined as the amplitude of the coherent peak divided by the coherent amplitude plus the twice thermal 
background at the peak center. This is in analogy with Michelson’s definition of optical fringe contrast. We define 
the final visibility as the average of visibility of the peaks in both the x- and y-directions. 

Reliable fitting in such a high dimensional parameter space is robust when the fitting routine starts with accurate 
initial guesses. To generate accurate initial conditions, the image is first low pass filtered to remove spatial frequencies 
corresponding to the diffraction peaks and higher, and then fitted to a single gaussian function. The coefficients 
resulting from this fit become the initial conditions of the thermal fraction in the final fit. Next, the guess for 
the thermal component is subtracted from the raw lineout and all remaining local maxima are identified with a 
peak finding algorithm. The peaks are sorted by amplitude, and the highest five locations are chosen to represent 
the initial position and amplitude in the initial conditions for the superfluid diffraction peaks. With these input 
conditions, nonlinear least-squares fitting of the six gaussian peaks gave reproducible fits with small residuals. An 
example of a typical fit is shown in Figure 1. 




FIG. 1. An example fit for 0,112 flux superfluid diffraction image after a hold time of 34 ms. The thermal fraction fit is shown 
in yellow with the components of the superfluid diffraction fit in cyan, purple, and green. Their sum is shown in dark blue laid 
over the raw data in orange. The residuals are shown in dark red. A visibility olV — 0.40 is derived from the amplitudes of 
this fit. 


Bandmapping Procedure 


To quantify the population imbalance in the two minima of the HH spectrum, we developed a band mapping 
technique, in which laser beams are not switched off instantly but instead ramped down to zero smoothly before 
time-of-flight absorption imaging. The requirements for band mapping in the a = 1/2 HH lattice are different from 
the standard bandmapping technique in two ways. First, the two Hofstadter sub-bands are connected at Dirac points 
and the total band gap is therefore zero, forbidding a globally adiabatic band mapping procedure. Since the atoms 
are condensed, they are sufficiently localized to quasimomentum states that are far away from the Dirac points and 
are locally gapped. Second, there is a strong potential gradient present, leading to Bloch oscillations in the lab 
frame. For low lattice depths, the Bloch oscillations would cause Landau-Zener transitions at the zone edge. This 
effect is minimized by ramping down the lattice relatively quickly, in about one Bloch oscillation time. In the case of 
topologically trivial superfluids, all the lattices are linearly ramped down in 0.43 ms. In the case of HH superfluids, 
Raman beams are linearly ramped down in 0.88 ms first to avoid excitations associated with the Raman drive, and 
then the lattices are linearly ramped down in 0.43 ms. 


Superfluid Diffraction Close to the Suerpfluid to Mott Insulator Transition 


As mentioned in the main text, we observe weak diffraction peaks in three-dimensional lattices for z-lattice depths 
up to 20 Er with X and y lattices held fixed at 11 Er^ demonstrating Bose-Einstein condensation in the HH lattice 
with strong interactions. The images in Figure 2 represent a selection of shots immediately after the ^35 ms ramp 
out of the Mott insulator state into the 1/2 flux superfluid. The doubling structure is clearly visible at lattice depths 
close to the Mott insulating transition for filling factors n> 2. 

No systematic studies were done for z-lattices greater than 11 E^ since much larger shot-to-shot fluctuations were 
encountered, probably due to larger sensitivity to technical noise sources. 
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FIG. 2. 1/2 flux superfluid diffraction pattern for increasing interaction strnength, controlled by the depth of the lattice in the 
z-direction. 


Collisional instability of the ground state in tilted and driven lattices 


The observed decay of a Bose-Einstein condensate in tilted and driven lattices is much faster than the decay 
in standard lattices. Although our lifetime seems to be limited largely by technical noise, we wish to discuss the 
limit given by the time dependent nature of the single-particle wave function representing the ground state of the 
Floquet Hamiltonian mm- Choudhury et al. study the stability of a condensate in a tilted one-dimensional lattice 
with resonant amplitude modulation using a high frequency perturbation analysis. The ground state of bosons in the 
corresponding Floquet Hamiltonian is a superfluid in a geometry with a lattice in one direction and free particle motion 
in the other two. The description of this system is simple yet nontrivial due to the metastability of the Wannier- 
Stark ladder and the time-modulated many-particle Floquet Hamiltonian. As pointed out in [1], one consequence 
of these two facts is that elastic collisions between atoms will deplete the condensate by using one or two quanta 
of the micromotion (equal to the offset energy A between neighboring sites) to create excitations. This will cause 
decoherence and relaxation to higher temperatures. 

There are two main processes contributing to this instability of the condensate. Wannier-Stark states - which are 
mostly localized in a given well n - have small components in lower wells n — 1, with amplitude oc 2J/A. This 
small component, which has energy A higher than the Wannier-Stark state localized in the well n — 1, can collide 
with atoms in the lower lattice site, thus creating two excitations with energy A/2 in the transverse direction (the 
creation of two exitations is necessary for conservation of transverse momentum) . The second mechanism involves 
the counter-rotating term of the phonon or photon assisted tunneling process. The co-rotating term enables resonant 
tunneling to a lower lattice site by stimulated emission of a phonon. In the counter-rotating process, the atoms tunnels 
to a lower lattice site, but absorbs a lattice phonon. Therefore, it has now an excess energy of 2A which can create 
two transverse excitations where each of them has energy A. 

These two-body elastic collisions convert modulation or tilt energy into motion along the transverse directions. The 
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characteristic relaxation times of the two processes are given by ri and T 2 '. 

^ k 


1 K‘^ 

T2 A"' 


^27r5(i;f-A). 


We can associate ri with a process that occurs even in the absence of modulation - the rate is proportional to the 
tunnel coupling squared between neighboring sites, whereas the rate of the T 2 process is proportional to the square of 
the phonon induced tunneling rate. Both decay processes are rather general and should apply to both tilted lattices 
and lattices modulated by superlattices. 

As a test of these predictions, we implemented this Hamiltonian experimentally using the same experimental 
procedure described in the main text for resonant amplitude modulation. We measured the coherence and lower band 
population lifetimes as shown in Figure 3. Approximate values in our experiment were = h x 550 — 830 Hz, = 
5.03 nm, d = 254 nm, K = h x 10 Hz, J = h x 30 Hz, and A = h x 3420 Hz, with gn as the chemical potential, 
the s-wave scattering length, K laser-assisted tunnelling rate, J the bare tunnelling rate, and A the energy tilt per 
lattice site. d~^ = f dx is the size of the Wannier statein the x-direction (j)x{x). These give us the values of 

l/r 2 = 0.003 — 0.005 s“^ and 1/ri = 0.13 — 0.20 s“^, which are longer than our observed lifetime.Uncertainty is driven 
by uncertainty in our density measurements due to redistribution during lattice ramp-up. However, during the long 
coherence lifetime, significant population loss was observed (due to inelatic two- or three-body collisions, or transfer 
to higher bands), and detailed further study is needed to identify the contributions of different processes to the decay 
rate. 


800 


s_ 

< 


E 


400 


—, ■ ' .I ' ' . I 

: T 

' .1 ' ■ .1 


o 

o 

o 

o 

^ 0 

i ® 

□ 

0 

o 


□ 

■ 

Z-1 _._...1_._...1_ 

□ i : 

I e: 

_...1_._...1_ 


0.3 


J' 

> 


30 

Hold Time (ms) 


300 


3000 


FIG. 3. Lifetime of a Bose-Einstein condensate in a ID lattice with resonant amplitude modulation. Shown are the decay of 
the number of trapped atoms, and the decay of the visibility of the diffraction pattern observed after ballistic expansion. 


The mechanism proposed by Choudhury et al appears to give a reasonable upper limit for collisional decay in 
ID lattices. We expect the same mechanism to explain at least part of the decoherence observed for both amplitude 
modulation in a two-dimensional lattice and the HH Hamiltonian. For a quantitative comparison, we modify the 
calculation in [1] to reflect that a 2D lattice has only one orthogonal direction of free motion. For resonant modulation, 
in the Floquet basis, the Hamiltonian is 

k ki,k2,k3 


with k 4 = ki + k 2 — ka, V as the system volume, and 5k and as the annihiliation and creation operators for states 
with quasimomentum k. The interaction parameter is given by 


9 = 
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where (t)x{x) and (t)y{y) are the Wannier functions in the tilted (x) and untilted (y) lattice directions, A is the lattice 
spacing, and d is the size of the Wannier state. The instantaneous single particle dispersion e\^{t) is given by 

^k{t) = —‘2K cos{kx) — 2Jy cos{ky) — 2K cos{kx — 2At) — 2Jx cos{kx — At) + ^ ^ . 

where Jx and Jy are the bare tunnelling rates in the x and y directions, and K is the effective tunnelling in the tilt 
direction. Following a similar argument, we find that the scattering rate is given by 

1 _ 1 1 

r T2 Ti 


1 

T2 


2{gn/2f 
Nh A2 


5]27r(5(4‘’)-A) 


1 

Tl 


2{gn/2f 4J2 

Nh A2 


^27r(5(4°)-A/2) 


with the effective dispersion: 


= 2K[1 — cos(/Ca,)] + 2J^[1 — cos{ky)] + 




2 


Note that ref. [T] uses a Hartree-Fock excitation energy with free particle dispersion. The high mean field energy in 
the one dimensional tubes created by a two dimensional lattice, gn^ is similar in energy scale to the tilt frequency, 
and we thus use the full Bogoliubov spectrum for the energy of quasi-particles. 

Since K and Jy are small, we neglect the dependence of on kx and ky^ giving 

p{u)=Y,2n5iEl^^-u) (2) 

k 



and thus 


1 _ gn Alias K‘^ A y 2 md‘^ 

T2 h d a \/ {gn)^ + ^^{gn)‘^ + 


1 _ gnAnas 4J^ A/2 y 2md‘^ 

Tl h d A ^ {gn)‘^ E (A/2)^ (^n) + \/ {gTi)‘^ + (A/^^ 


For our system, representative numbers are gn = h x 740 — 1740 Hz, and A = h x 3420 Hz, = h x 10 Hz, 
J = h X 30 Hz, as = 5.03 nm, d = 254 nm, and m = 1.4 x 10“^^ kg, resulting in 1/t2 = 0.006 — 0.014 s“^ and 
1/ri = 0.29 — 0.67 s“^. This is much smaller than the decay rates measured in the experiment. 
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FIG. 4. Lifetime of the total atom number (red) in the lowest band of the 3D lattice as a function of the z-lattice depth. For 
comparison, the coherence lifetime (blue) as a function of lattice depth is shown on the same plot. In the main section of the 
paper the coherence data is displayed as a function of the calculated Hubbard interaction parameter, U. 


Loss Processes in the 3D HH Model 

By monitoring the loss of atoms we can show that certain processes are negligible for the loss rate of the superfluid 
state. Figure 4 shows the lifetime of the population in the lowest band of our lattice - or equivalently the total atom 
number lifetime, since the strong tilt will quickly remove population promoted to higher bands due to Landau-Zener 
tunneling to the continuum. The lifetime for atom loss from the lower band is about an order-of-magnitude longer 
than the coherence lifetime and shows a slight decrease with increasing lattice depth. For comparison, a BEG in an 
untilted lattice exhibits coherence times exceeding one second in deep lattices. 

We conclude that the population loss seen in Figure 4 is dominated by some process beyond heating in standard 
lattices (spinflip collisions, magnetic field noise, lattice intensity noise, and beam pointing noise). However, the 
dominant decay process for the BEG in the HH lattice are excitations within the lowest band. 

A 3D lattice should be helpful for reducing collisional heating. Qualitatively, in a Eermi’s golden rule picture of 
heating, the decay rate depends on the dispersion relation of the transverse modes [1]. Thus, if a 3D lattice has a 
bandgap at the energies A and A/2, at which quasiparticles can be created, we would expect to see a “gapping out” of 
collisional heating and a corresponding increase in the lifetime upper limit. Since we do not see a qualitative change in 
lifetime when changing the bandgap, we conclude that, similarly to a 2D lattice, technical fluctuations are the likely 
limit of the lifetime. 
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